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ABSTRACT 


Prediction  Intervals  for  future  observations  in  serially 
correlated  samples  from  a  normal  distribution  are  derived. 
The  results  are  extended  to  predict  a  future  observation  in 
a  linear  trend  model.  The  properties  of  the  prediction 
intervals  are  examined. 


TABLE  OF  CONTENTS 


I.  INTRODUCTION  .  5 

II.  PREDICTION  INTERVALS  FOR  NORMAL  DISTRIBUTION  .  7 

A.  SUMMARY  OF  KNOWN  RESULTS  .  7 

B.  DEFINITIONS,  NOTATIONS  AND  REVIEW  OF  SIMPLE 

LINEAR  REGRESSION  ANALYSIS  .  8 

C.  STATIONARY  GAUSSIAN  MARKOV  PROCESSES  .  10 

III.  PREDICTION  INTERVAL  FOR  A  FUTURE  OBSERVATION 

IN  SERIALLY  CORRELATED  NORMAL  SAMPLES  .  12 

A.  DERIVATION  OF  THE  PREDICTION  INTERVAL  .  12 

B.  PROPERTIES  OF  THE  PREDICTION  INTERVAL  .  16 

C.  DISTRIBUTION  OF  *F* .  20 

D.  NUMERICAL  EXAMPLE  .  22 

IV.  SIMULATION  OF  SERIALLY  CORRELATED  SAMPLES  .  26 

A.  SIMULATION  FROM  STANDARD  NORMAL  SAMPLES  .  26 

B.  TESTING  THE  SAMPLES  FOR  NORMALITY 

AND  CORRELATION  .  30 

V.  RESULTS  OF  SIMULATION  .  32 

VI.  APPLICATION  TO  LINEAR  TREND  MODEL  .  35 

A.  PREDICTION  INTERVAL  FOR  X2k+1  .  36 

B.  NUMERICAL  EXAMPLE  .  40 

COMPUTER  OUTPUT  .  43 

PLOTS  OF  INTERVALS  .  58 

COMPUTER  PROGRAMS  .  63 

LIST  OF  REFERENCES  .  69 

INITIAL  DISTRIBUTION  LIST  .  70 


4 


I .  INTRODUCTION 


\ 

A  prediction  interval  is  an  interval  that  contains  a 
future  observation  with  a  pre- specified  probability.  The 
limits  of  the  interval  are  functions  of  known  observations 
from  a  family  of  known  distributions  with  given  parameters. 
Though  prediction  intervals  resemble  confidence  intervals , 
they  differ  from  the  latter  conceptually.  A  confidence 
interval  covers  the  value  of  a  parameter  of  a  distribution. 

A  prediction  interval,  on  the  other  hand,  encloses  the  value 
of  a  random  variable. 

A  prediction  interval  for  a  single  future  observation 
and  simultaneous  intervals  for  independent  and  identically 
distributed  samples  from  N(y,cr2)  were  derived  by  Hahn  [2]. 

For  correlated  samples  with  a  prescribed  correlation  struc¬ 
ture,  prediction  intervals  were  developed  by  Choi  [3].  It 
was  shown  by  him  that  many  of  Hahn's  prediction  intervals 
are  valid  even  for  the  case  where  the  samples  were  correlated 
with  the  specified  correlation  structure. 

In  this  thesis,  the  observations  are  assumed  to  belong  to 
a  multivariate  normal  family  with  mean  vector 


y  =  (y  ,y,  .  .  .  ,y)T 
nxl 

and  correlation  matrix 


-  5  - 


(1-1) 


Let  Xt(t=. . . -2, -1,0, 1, 2. . . )  be  a  stationary  Gaussian 

markov  process  of  order  one  with  E(Xt)  =  y  for  all  t  and 

2  k 

the  covariance  between  Xfc  and  =  a  p  for  all  t  and  k. 

Stochastic  processes  of  this  type  are  used  to  model  certain 
types  of  oceanographic  and  metereological  phenomena  and 
some  aspects  of  the  stock  market.  Any  finite  sub  process 
will  then  belong  to  the  above  multivariate  normal  family. 

Chapter  2  deals  with  prediction  intervals  from  a  normal 
distribution  and  contains  a  summary  of  known  results  and 
a  review  of  general  linear  hypothesis.  In  chapter  3  a 
prediction  interval  for  a  future  observation  in  a  first  order 
markov  process  is  derived  and  the  properties  of  this  interval 
are  examined.  The  effect  of  various  parameters  such  as 
variance,  correlation  and  sample  size  on  the  prediction 
interval  are  also  discussed  in  this  chapter.  Chapter  4 
contains  a  discussion  of  a  simulation  to  generate  serially 
correlated  random  variables,  which  are  used  to  empirically 
verify  the  theoretical  conclusions  drawn.  The  results  of 
the  simulation  are  presented  in  chapter  5.  Extensions  of 
the  results  to  a  linear  trend  model  are  discussed  in  chapter  6. 
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II.  PREDICTION  INTERVALS  FOR  NORMAL  DISTRIBUTION 


A.  SUMMARY  OF  KNOWN  RESULTS 

Consider  a  random  sample  of  size  N=n^+n^  drawn  from  a 

2 

normal  distribution  N(y,o  ).  The  samples  n^  form  the  group 

of  initial  samples  and  n^  form  the  group,  called  future 

_  2 

samples.  Based  on  the  mean  X  and  variance  S  of  the 

ni  ni  2 

initial  samples,  prediction  intervals  for  X^  and  S^  were 
obtained  by  Hickman  [6].  Hahn  [2]  derived  prediction 
intervals  for  the  mean  and  variance  of  the  second  sample 
and  also  simultaneous  prediction  intervals  for  the  variance 
of  each  of  k  additional  random  samples  of  size  n^  based  on 
the  information  obtained  from  the  initial  sample.  Hahn's 
[2]  two  sided  100  r  %  prediction  interval  to  contain  a 
single  additional  observation  X  is  given  by 

Y  ±  t[n-l;(l+r)/2]  (l+^)%  •  S 

where  Y  is  the  sample  mean,  S  is  the  sample  variance  and 
t[n-l , (1+r) / 2 ]  is  obtained  from  the  tables  of  t  distribution 
with  (n-1)  degrees  of  freedom.  A  100  r  %  simultaneous 

prediction  interval  to  contain  X^,  .  Xk>  k  future 

observations  is  given  by 

Y  +  R(k,n,r)S 

where 

R(k,n,r)  =  U(l+^)% 

The  value  of  the  function  U  can  be  obtained  from  the  special 
tables  given  in  reference  2. 
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A  prediction  interval  for  a  single  future  observation  as 
well  as  simultaneous  intervals  for  a  specified  number  of 
future  observations  have  been  derived  by  Choi  [3]  for  the 
case  where  the  observations  are  correlated  and  belong  to  a 
multivariate  normal  distribution  with  mean 

T 

y  =  (y  ,y . . y)  , 

and  covariance  matrix , 

V  =  %(  H  +  H1)  +  ct(  I  -  E  ) 
nxn  nxn  nxn  nxn  nxn 


h^(i=l,2 . n)  and  a  are  positive  constants.  E  is  an  nxn 

matrix  all  of  whose  elements  are  unity. 

B.  DEFINITIONS,  NOTATIONS  AND  REVIEW  OF  SIMPLE  LINEAR 
REGRESSION  ANALYSIS 

A  brief  review  of  the  results  and  notation  in  a  simple 
linear  regression  model  are  presented  below.  These  results 
are  used  in  chapter  3  in  deriving  prediction  intervals  for 
a  first  order  markov  process.  By  way  of  notation,  a 
capital  letter  with  an  underbar  represents  a  vector  or  a 
matrix.  Consider  the  simple  linear  regression  model 

Y  =  X  B  +  e 
nxl  nx2  2x1  nxl 

2 

where  e  v  N(0,a  1^) 


(2-1) 


which  in  scalar  notation  can  be  written  as 


Yi  =  a  +  bXi  +  ei  i  -  1,2, . . . ,n. 

The  least  squares  estimator  for  B  is  obtained  by 

Min  Q  =  (Y  -  XB)T(Y  -  XB ) , 

B 

that  is  3  =  S_1XTY  , 

where  S  =  XTX  , 

and  B  *  N(B,  a2S_1) .  (2- 

T 

B  is  the  column  vector  (a,b) 

It  is  well  known  that  the  random  variable 

Y  ~  N (XB ,  a2I) 

EX2  -EXi 

i  f 

EXi  n 

where 

n  -  2 
SXX  =  E  (X.  -  X)^  . 

i=l  1 

^2  2 
Let  a  be  the  usual  unbiased  estimator  of  a  ,  given  by 

a2  =  (SYY  -  b2SXX) / (n  -  2)  ,  (2- 

n  -  2 

where  SYY  =  Z  (Y,  -  T)  . 
i-1  1 


9 


Let  Y  be  the  predicted  value  of  the  random  variable  at  a 
specified  value  of  the  independent  variable  x. 

A  100(l-a)7o  prediction  interval  for  the  dependent  vari¬ 
able  Y  can  be  obtained  from  the  relation 


r 

A 

Y-Y 

<  t 

=  1  -  a 


(2-4) 


Equation  (2-4)  will  be  used  in  deriving  a  prediction  interval 
for  x  in  chapter  3 . 

C.  STATIONARY  GAUSSIAN  MARKOV  PROCESSES 

Let  Xt  (t  =  . . . -2 , -1 , 0 , 1 , 2 , . . . )  be  a  stationary  gaussian 
markov  process,  i.e.,  an  autoregressive  stochastic  process, 
of  order  one  with  E(Xt)  =  m  and  covariance  between  Xt  and 


X 


t+k 


2  k 
a  p  . 

To  derive  the  prediction  interval  for  a  single  future 


observation  based  on  the  observations  up  to  the 

results  of  Ogawara  [1]  will  be  used.  He  has  shown  that 
when  the  values  of  X2k-1  (k=l >  2  ,  •  •  ,n+l)  are  fixed, 

(k=l , 2 , . . . ,n)  are  mutually  independent.  The  conditional 
probability  densities  are  given  by 


f(X2k  X2k-1’  X2k+1) 


EXP 


no 


_1 

2a 


where 


a  =  m(l-p2) / (1+p2) 

b  =  2p/(l+p2) 


x2k-  (a+bxrPi 


(2-5) 
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A  prediction  interval  for  X2^+^  will  be  derived  in  chapter 
3  based  on  this  conditional  distribution  of  the  even  numbered 
observations  conditioned  on  the  odd  numbered  ones. 


III.  PREDICTION  INTERVAL  FOR  A  FUTURE 
OBSERVATION  IN  SERIALLY  CORRELATED  NORMAL  SAMPLES 


A  prediction  interval  to  contain  a  single  future  obser¬ 
vation  for  samples  from  a  multivariate  normal  distribution 
with  a  correlation  matrix  of  the  form  defined  in  equation 
(1-1),  is  derived  in  this  chapter.  In  section  A  a  conditional 
prediction  interval  to  contain  a  single  additional  observation 
based  on  the  conditional  distribution  of  X2^  is  obtained. 
Section  B  deals  with  the  properties  of  the  prediction  interval. 
In  section  C  the  dependence  of  the  prediction  interval  on 
sample  size,  p  and  sigma  are  discussed. 


A.  DERIVATION  OF  THE  PREDICTION  INTERVAL 

Assume  a  sample  of  observations,  X^,  X^,  X^ . X-2^ 

X2k  from  a  normal  distribution  with  mean  zero  and  covariance 
matrix  V  such  that  Cov  (,Xt>  xt+^)  =  P  • 

The  conditional  probability  densities,  as  indicated  in 
equation  (2-5),  are 

-  iHlx2k-  <a+bx;>l 2 

O  ^  J 

with  conditional  mean 


f(X 


2k 


X2k-1’  X2k+1) 


\/  2ira 


EXP 


E (X2k  Xk)  =  a+bXk  ’ 


and 


K  =  (X2k-l  +  X2k— L)/2 


(3-1) 


E  (X 


2k 


*k> 


considered  as  a  function  of  X^  is  callled  the 


regression  function  of  X2^,  on  X£.  Graphically,  it  represents 


the  locus  of  the  center  of  gravity  of  the  conditional  random 


variables  X2k J  as  a  function  of  X£. 

As  a  result  of  Ogawara’s  Theorem  the  following  conditional 
simple  linear  regression  model  may  be  assumed, 


X2k  a  +  bXk  +  ek  ’  k  1,2,  .  .  . ,  n 


(3-2) 


where  =  N(0,  a  ). 

The  maximum  likelihood  estimators  of  the  parameters  a,  b 
and  a ^  are  given  by 

b  =  SXY/SXX 


a  ’  x2k  '  b  *k 

a2  -  (SYY  -  b2  SXX)/(n-2)  , 


where 


=  1  /  xl^x3  +  X3^~X5  +  +  x2k-l*x2k+l 

k  n  \  2  2  '  ‘  '  2 


n  \  x2  +  x4  +  '  "  +  x2k 


k=i  ( X2k  f  'n  ^ 


k£1  K  X2k  •  n  xk  x2k 


^(^k)2  -n  x2k2  • 


(3-3) 


At  a  known  value  of  x£  the  predictor  x2k  is  given  by 


x„,  =  a  +  b  x,' 


(3-4) 


A  100(l-ct)7o  prediction  interval  for  the  dependent  variable 
x2k  can  obtained  from  the  relation 


r , 


From  the  above  relation  a  prediction  interval  for  x2k+1  can 
be  obtained  noting  that  x£  is  a  linear  function  of  x2k+1 . 
The  inequality  inside  the  brackets  can  be  written  as 


(3-6) 


Substituting  for  x2k  -  (x2k  -  b  x£)  +  b  x£  and  squaring  both 
sides,  the  inequality  (3-6)  reduces  to 


2b(x2k-  x2k)(xk  '  Xk} 


t2o2(l  +  i)  +  t2a2(xk  - 


xk>  2 


xo'/sxx 


or 

(fa2  *  TnBr)(*k  ‘  )2  ■  2^(x2k  ■  x2k)(xk  *  xk) 

+  u2k-x2k>2  -  t2S2a  +  i)<  o  (3-7) 

This  is  a  quadratic  equation  of  the  form 

ACxk  ‘  xk)2  +  B(xk  -  *£)  +  C  <  0  .  (3-8) 

14 


i 


where 


A  = 


.2~2 
t  a 

sxx 


B  =  -2b  (x2k  -  x2k) 


C  (x2k  "  *2k)2 


t2a2  5±i 


(3-9) 


Thus,  the  probability  statement  in  equation  (3-5)  is  equiva¬ 
lent  to 


P[A(xk  "  x£)2  +  B(x^  -  x£)  +  C  <  0]  »  1  -  a 


(3-10) 


A  prediction  interval  for  xk  can  now  be  obtained  in  terms  of 
the  roots  of  the  quadratic  equation  (3*8)  viz., 


(K  -  - 


-B±/ B2  - 


4AC 


2A 


the  term  under  the  radical  sign  simplifies  to 


2  9  ~  9 

Bz  -  4AC  =  4tZa 


-  x  '2 
"2k  x2k/ 


SXX 


,  n+1  f C2  t2a2 

+  —  Vb  -  W 


and  hence 


-  =.  .  £(x2k'x2k 
Xk  *k  +  - 


(x2k'x2k)2  +  n+1 


SXX 


'b2 


\ 


n 


b2-  C  a 


Since  x^  =  (x2k_1  +  x2k+l)/2  * 


(3-11  ) 


Sxx  /  (3-12) 


The  unknown  future  value  x2^+l  can  be  written  in  terms  of  the 
known  value  as 

D  ±  E  where 


2k+l 


15 


(3-13) 


C2  t2o2 


D  “  2xk  "  x2k-l  +  2b  (x2k  "  x2k^b  “  ~§XX"  ) 


E  =  2ta  \/  (x2k  -  x2k)2/SXX  +  ^(b2  b 


t2a2x  /f  2  t2g2\ 

sxx '  /  \D  “  ~sxT / 


(3-14) 


Thus  a  prediction  interval  for  x2k+-^  is  given  by 


D  -  E  <  x2k+]_  <  D  +  E  , 


(3-15) 


where  D  and  E  are  given  by  equations  (3 -13)  and  (3-14)  respec¬ 
tively.  The  properties  of  this  interval  are  discussed  in 
section  B. 


B.  PROPERTIES  OF  THE  PREDICTION  INTERVAL 

The  prediction  interval  in  the  previous  section  was 
obtained  by  solving 

A(xk  -  +  B(xk  -  +  C  =  0 


where  the  coefficients  A,  B  and  C  are  defined  in  (3-9).  The 
roots  of  the  equation  and  hence  the  form  of  the  prediction 
interval  will  depend  on  the  relative  sizes  of  A,  B  and  C. 

The  effect  of  these  coefficients  on  the  prediction 
interval  will  now  be  examined  breaking  xn  to  four  ail  in¬ 
clusive  cases. 


Case  1. 


2 

A  >  0  and  B 

"2 

A  >  0  b 


-  4AC  >  0 

+2*  2 
t  a 

>  “sxx 


16 


or  b2SXX 
t  o 


>  1  . 


/VO 

b^SXX 
t  a 


(3-16) 


Bz  >  4AC  implies  that 
4  b2  Cx2k  -  x2k)2  >  4  (b2 


4  (b^  - 


i-K  2 

to  _  2 

Tx3T)  <x2k  '  x2k) 

t2o  2  t2o  2 

-s»r)  — —  (n+1) 


or  \  b 


t2a  2\  t2;  2 


2 

t  a 


"5  (n+1)  >  §xx~  <x2k  “  x2k^ 


"2 

b  SXX 


-  1  > 


n  (x2k  '  x2k)2 


t  a 


-2 

b^SXX 

2"2 
t  a 


>  1  - 


n  (x2k  '  x2k)2 


t  a 


To  satisfy  the  conditions  A  >  0  and  a  >  4AC  it  is  necessary 


that 


„  ,  ,  „  .  ,  n  <*2k  -  x2k>2 

F  >  1  and  F  >  1  -  - 2 — 


(3 -16a) 


—  2  2  ~  2 

Since  the  quantity  (x2k  -  x2k)  /t  a  is  always  positive  the 
conditions  in  equations  (3 -16a)  are  equivalent  to  the  condition 
F  >  1.  In  this  case  the  prediction  interval  for  x2k+1  will 


17 


be  of  the  form 


[(D  -  E)  ,  (D  +  E)]  , 


This  type  of  two-sided  interval  is  referred  to  as  TYPE  1, 
where  D  and  E  are  as  defined  in  (3-13)  and  (3-14)  respectively 


Case  2.  Is  -  4AC  >  0  and  A  <  0 


A  <  0 


.  2 
t  a 

SXX 


<  0 


or  b  < 


c2; 2 

SXX 


or 


b2SXX 
t  a 


i.e.,  A<0  F  <  1, 


(3 • 17) 


„2  n  (x2k  '  x2k) 

B  -  4AC  2.  0  -  F  - -n - 

t  a 


(3-18) 


Thus ,  if  F  <  1  and  F  2.  1  - 


n  (x2k  "  x2k)2 


t  a 


the  resulting  interval  will  be  of  the  form, 

<-  «,  D  +  E]  [D  -  E,  «  )  (3-19) 

and  will  be  called  a  TYPE  2  interval. 

2 

Case  3.  If  B  -  4AC  <  0  and  A  ^  0  the  interval  is  called 


TYPE  3. 


The  roots  of  the  quadratic  are  complex  and  the  prediction 

interval  must  be  taken  as  (-  00 ,  +  °°)  .  This  situation  occurs 
2 

when  B  -  4AC  <  0  . 


CASE  4.  A  =  0 

A  -  0  *  1  .  (3.21) 

t  a 

Thus,  when  F  =  1,  the  interval  is  of  type  4  and  is  of  the 
form  [G,®  )  , 

where  G  =  2x£  -  x2k_1  -  b  SXX  (**±1)  +  X^  a  X2k  ,  (3-22) 

To  summarize  the  results,  the  quadratic  equation  that 
leads  to  the  required  prediction  interval  is  of  the  form 

A(x£  -  x£)2  +  B(x^  -  ££)  +  C  -  0, 

and  four  different  types  of  intervals  can  result  depending 
on  the  val  -r»  of 

„  _  b  SXX 

F  "  .  2A2  • 

t  a 


The  four  cases  are 


i)  when  F  >  1,  the  interval  is  of  the  form 


ii) 


iii) 


iv) 


[(D-E) . (D+E)3 - 
when  F  <  1  and  F  z.  1- 


n  (x2k  '  x2k)2 


n+1 


JT 

t  a 


the  interval 


is  of  type  2  and  is  of  the  form  (-  «>,  D+E] ,  [D-E ,°°) 
(x. 


when  F  <  1- 


n 


L2k 


H+r 


-  x2k>2 


TTT 
t  a 


the  interval  is  of 


type  3  and  is  of  the  form  (-oo^+oo) 

when  F  =  1,  the  interval  is  of  type  4  i.e.,  [G,“). 


At  a  given  level  of  significance  a,  from  the  relation 
(3-5)  and  the  conclusions  drawn  in  this  section,  the  following 
relation  should  hold: 


4 

E  Pr  [occurrence  of  type  i  interval] 

=1  th 

x  Pr  [i  interval  contains  the  predicted  value] 

=  1  -  a  (3-23) 


The  above  conclusions  are  verified  with  simulated 
samples  in  chapter  5. 

Section  C  examines  the  distribution  of  the  quantity  F 
defined  in  equation  (3-16). 

C.  DISTRIBUTION  OF  ’F’ 


The  occurrence  of  any  particular  type  of  interval  depends 
on  the  value 


F  = 


b2  SXX 

2A2 

t  a 


The  distribution  of  F  can  be  shown  to  be  a  noncentral  F 
distribution  as  follows. 
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I 


*  9 

From  (2*2),  it  can  be  shown  that  b  is  N(b,a  /SXX) 

A 

or  9,  N(0, 1)  , 

o//Sl QC 


%  x2(D, 

az/SXX 


and  b2-^SXX  ^  x2  (1>  x)  _  (3-24) 

a 


where  X 


b2  SXX 

2 

a 


is  the  noncentrality  parameter. 
From  (2-3)  the  distribution  of 


(n-2)oi 


is  X  (n-2) 


(3-25) 


The  ratio  of  (3-24;  and  (3-25)  devided  by  their  respective 
degress  of  freedom  results  in  a  noncentral  F  [5]  statistic 
with  degrees  of  freedom  (1,  n-2)  and  noncentrality  parameter 
X .  That  is 


b  •  SXX  J2^  .  ^  _  0 

- xtj —  =  t  F  v  IF  (1,  n-2,  X)  , 


and  the  noncentrality  parameter 
,  b2  SXX 


(3-26) 


Substituting  for  b  from  equation  (2-5),  the  noncentrality 
parameter  X  can  be  written  as 
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n  o 


4P 

(1+P2)2 


sxx 

~2 

a 


(3-27) 


D.  NUMERICAL  EXAMPLE 

The  Model  developed  in  section  A  is  tested  on  the  data 
obtained  from  Dow-Jones  monthly  averages  to  predict  a  future 
monthly  average,  based  on  the  averages  of  the  previous  months. 
In  most  of  the  cases,  the  future  value  was  contained  in  the 
prediction  interval.  To  illustrate  an  example  the  data  for 
the  years  1966  and  67  are  given  below.  The  prediction 
intervals  are  computed  with  the  help  of  tne  APL  program 
PREDICT,  whose  listing  is  given  on  page  68. 


1966  1967 


Month  Ended 

D  -  J  Ind 

Month  Ended 

D  -  J  Ind 

January 

31 

983.51 

January 

31 

849.87 

February 

28 

951.89 

February 

28 

839.37 

March 

31 

924.77 

March 

31 

865.98 

April 

29 

933.68 

April 

28 

897.05 

May 

31 

884.07 

May 

31 

852.56 

June 

30 

870.10 

June 

30 

860.26 

July 

29 

847.38 

July 

31 

904.24 

Augus  t 

31 

788.41 

August 

31 

901.29 

Sept 

30 

774.22 

Sept 

29 

926.66 

October 

31 

807.07 

October 

31 

879.74 

November 

30 

791.59 

November 

30 

875.81 

December 

30 

785.69 

December 

29 

905.11 
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The  above  data  are  stored  in  the  variable  DMA.  The  pro¬ 
gram  Predict  is  run  with  different  sample  sizes  and  the  out¬ 
put  is  shown  on  pages  24  and  25. 

Using  14  monthly  averages  from  January  '66  to  February 
'67,  prediction  intervals  tor  March  '67  are  computed  and  the 
exact  value  of  865.98  is  contained  in  the  interval. 

Using  given  values  up  to  April  '67,  the  prediction  in¬ 
tervals  for  May  ’67,  also  contain  the  true  value.  With  a 
sample  of  size  18,  the  prediction  intervals  for  the  nine¬ 
teenth  month  i.e.,  July  '67  are  computed  and  they  contain  the 
true  value.  Though  the  average  of  September  '67  falls  outside 
the  interval,  the  average  of  November  ’67  is  contained  within 
the  interval. 

The  reduction  in  the  length  of  prediction  interval  with 
increase  in  sample  size,  is  supported  by  the  above  example. 

The  data  are  tested  for  the  required  correlation  structure 
and  found  to  fulfill  the  requirements  partially.  The  statistics 
of  the  data  are 

Mean  =  870.83 

Standard  deviation  =  55.0 

Corr.  Coefficient  p  =  0.82 


PREDICT  14+ DMA 


THE  F  STATISTIC  IS 
THE  LOWER  PREDICTION  LIMIT 
THE  UPPER  PREDICTION  LIMIT 
THE  LENGTH  OF  THE  INTERVAL 
TRUE  VALUE  OF  THE  AVERAGE 


PREDICT  16+ DMA 


THE  F  STATISTIC  IS 
THE  LOWER  PREDICTION  LIMIT 
THE  UPPER  PREDICTION  LIMIT 
THE  LENGTH  OF  THE  INTERVAL 
TRUE  VALUE  OF  THE  AVERAGE 


PREDICT  18+DMA 


THE  F  STATISTIC  IS 
THE  LOWER  PREDICTION  LIMIT 
THE  UPPER  PREDICTION  LIMIT 
THE  LENGTH  OF  THE  INTERVAL 
TRUE  VALUE  OF  THE  AVERAGE 


PREDICT  20+ DMA 


THE  F  STATISTIC  IS 
1HE  LOWER  PREDICTION  LIMIT 
THE  UPPER  PREDICTION  LIMIT 
THE  LENGTH  OF  THE  INTERVAL 
TRUE  VALUE  OF  THE  AVERAGE 


4.132073877 

727.1128824 

1080,631201 

353.5183188 

865.98 


5.489038103 

598.3935849 

900.5564036 

302.1628187 

852.56 


5.1305x7543 

708.7039 

1009.101896 

300.3979963 

904.24 


5.777973921 

573.5929826 

864.8705316 

291.277549 

926.66 
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PREDICT  22+DMA 


THE  F  STATISTIC  IS 
THE  LOWER  PREDICTION  LIMIT 
THE  UPPER  PREDICTION  LIMIT 
THE  LENGTH  OF  THE  INTERVAL 
TRUE  VALUE  OF  THE  AVERAGE 


6.994239226 

633.1176654 

899.7280383 

266.6103729 

875.81 
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IV.  SIMULATION  OF  SERIALLY  CORRELATED  SAMPLES 


To  study  the  effect  of  p ,  a  and  n  on  the  type  of  pre¬ 
diction  interval  obtained  it  is  necessary  to  generate  random 
variables  with  the  required  variance  covariance  matrix. 

These  can  be  generated  as  linear  combinations  of  standard 
normal  random  variates. 

A.  SIMULATION  FROM  STANDARD  NORMAL  SAMPLES 

Suppose  Zp  Z^,  ■  • • ,  Zn  are  random  samples  from  N(0,  I). 
To  get  serially  correlated  random  variates  X1  ,  X9,...,Xri, 
let  X]_ 


The  coefficients  C-^ ,  . ...C^  are  then  determined  by  solving 
equations  of  the  form, 

Ir 

Cbv  (X  ,  X^^)  =  p  for  all  t  and  k  from  1  to  n. 

Consider  the  first  equation 


E(XX)  =  E(ZL)  =  0 

Var  (X1)=  Var  (Z-^  =  1  (4 ’2) 


C^Z^+C2Z2 > 


C’Zt+C1Z0+ . +C'Z 

112  2  n  n 


(4-1) 
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(4-3) 


From  the  second  equation 

e(x2;  =  E(c1z1+c2z2) 


=  C1E(Z1)+C2  E(Z2)  =  0 


since  Z^  and  Z2  are  independent; 
Var  (X2)  =  Var (C1Z1+  C2Z2) 


Var  Z^  4-  C2  var  Z2  j  i.e.  , 


Cl2  +  c22  *  1 


(4-4) 


Cov(X^ ,  X2;  =  E  [(X1  -  E(X1))(X2  -  E(X2))] 

=  E  [X1-  X2] 


=  E  [Z1- (C1Z1+C2Z2)] 

—  E  [C^Z-^  +C2Z1^2]  j  i.e, 
C1E(Z12)+C2E(Z1Z2)  =  p  . 

Since  Z1  v  N(0,  1)  ,  Zj2  v  X2(l). 


also  since  Z^  and  Z2  are  independent,  the  above  equation  can 
be  simplified  to 

=  p 


(4-5) 


From  4.4  by  substitution  we  get  C2  =  y  1  - 


or  Xr 


pzl  +  y  i  -  p1 


(4-6) 


or  X, 


PXX  +  -  p: 


(4-7) 
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Similarly  it  can  be  shown  that 


=  Px2  +  V  1  -  P: 


(4-8) 


The  result  can  be  generalised  in  the  form 

xi.  -  pxi-i +  V 1  - p;!  zi 

i  =  2,  . . . ,  n. 


(4-9) 


or  Xi  =  p1'1  Z]_  +  p1"2^  1  -  p2  Z2  + - +  -  p2  Zt  (4-10) 

i  =  1,  . . . ,  n. 

The  results  obtained  in  (4-10)  can  be  expressed  in  the  follow¬ 
ing  matrix  notation: 


X 

=  A  Z 

nxl 

nxn  nxl 

where 

X 

nxl 

=  cxlf  x2>  .... 

VT- 

Z 

nxl 

=  (Zx ,  Z2 ,  . . .  > 

Zn>T  • 

and  A 

is  a 

nxn  matrix  of  the 

form 

(4-11) 


A  = 


0 

/TT 7 


0 

0 


“-1  -"-2  /TT7 


A  -  p2 _ o 

. ATTT7 


(4-12) 


It  a  random  variable  Z  is  distributed  normally  with 

nxl 

mean  y  and  covariance  V  and  if  A  is  an  nxn  matrix  of  rank  n, 


then  the  vector 
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Y  =  A  Z 

nxl  nxn  nxl 

is  distributed  as  multivariate  normal  with  mean  vector  A  y 
and  covariance  A  V  AT.  As  a  result,  the  random  variable  X 
(4-11)  is  distributed  as  multivariate  normal  with  the  re¬ 
quired  covariance  matrix  V.  Here 

Z  v  N(0,  I)  , 

where  I  is  an  nxn  identity  matrix. 

Hence,  X  %  N(0,  A  I  AT)  . 

T  T 

We  now  show  that  A  1^  A  =  A  A  is  a  matrix  of  the  required 
form  (for  n  =  4) 

10  0  0 

p  /l-p2  0  0 

p2  p/l-p2  /l-p2  0 

3  2  l  1  l  I  Z  1 

p  p  /l-p  p/l-p  /l-p 

1  p  p2  p3 

o  /i/p2  p/i-p2  p2/i-p2 

0  0  /l-p2  p/l-p2 

ooo  /i-p2 


p3+p(l-p2) 

p4+p2(l-p2) 


2 

P 

p3+P (1-p2) 
p4+p2(l-p2)+(l-p2) 
p3+p3 (l-p2)+p (1-p2) 


4,  2/1  2N 

p  +p  (1-p  ) 
p5+p3(l-p2)+p(l-p2) 
p^+(l-p2) (p4+p2+l) 


which  reduces  to 

“  2  3 

1  P  P  P 

,  2 

p  1  p  P 

A  At  =  2 

P^  P  1  P 

_P3  p2  p  1 

The  above  result  can  be  generalized  for  any  value  of  n. 

B.  TESTING  FOR  NORMALITY  AND  CORRELATION 

The  simulated  samples  are  tested  for  normality  and  the 
correlation  structure.  The  simulated  samples  are  tested  by 
standard  Kolmogorov- Smirnov  test.  The  results  showed  that 
samples  do  appear  to  be  from  the  specified  normal  distribution. 
The  correlation  between  and  is  given  by 

Cov(X. ,  X . ) 

R  =  — -  -L  2  -  (4  •  13) 

War  X.  •  Var  X. 
i  J 

Confidence  limits  for  correlation  coefficient  p  are  developed 
using  the  approximation. 
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Z  =  %  In  Q  -  N<%  In  ^  ,  jjij) 

A  957o  confidence  interval  for  p  is  given  by 


(e2a-l)/(e2a+l)  <  P  <  U2b-l)/(e2b+l), 


(4-14 


where 


a  =  Z  -  and  b  =  Z  +  if— 

/n-3  /n-3 


The  simulated  samples  are  also  tested  for  the  required 
correlation  structure  by  examining  the  confidence  limits  for 
correlation  between  (X^,  X ^ )  using  equation  (.4-14).  The 
results  showed  that  the  samples  do  appear  to  have  the  pre¬ 


scribed  correlation  structure. 
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V.  RESULTS  OF  SIMULATION 


In  this  chapter,  using  simulated  samples,  the  probabili¬ 
ties  of  occurrence  of  various  types  of  intervals  and  their 
dependence  on  the  value  of  p ,  n  and  a  are  studied.  The  pro¬ 
gram  PROBS  FORTRAN,  whose  listing  is  given  on  pages  63-66 
computes  the  probabilities  for  the  four  types  of  intervals. 

For  given  values  of  p ,  n  and  o  the  program  generates  one 
thousand  samples  and  records  the  frequency  of  occurrence  of 
each  type  of  interval.  In  each  case  the  program  tests  whether 
the  predicted  value  is  contained  within  the  interval.  Finally, 
for  a  given  level  of  significance  a  the  relation  (3-23)  is 
verified. 

Choosing  a  level  of  significance  a  =  0.05,  the  simulation 
was  run  for  values  of  p  ranging  from  0.1  to  0,9,  for  sample 
sizes  from  10  to  50  and  for  sigma  varying  from  1  to  5 .  For 
a  =  0.1  and  0.02  the  program  is  executed  to  study  the  effect 
of  the  level  of  significance  on  the  results. 

The  computer  output  is  shown  on  pages  43  to  57.  Tables 
I  to  VI  show  the  variation  of  probabilities  of  types  of 
intervals  with  change  in  a  at  different  sample  sizes,  keeping 
a  and  a  constant.  Tables  VII  to  IX  show  the  variation  of 
probabilities  of  types  of  intervals  with  a  change  in  p  at 
values  of  a  =  0.1.  Tables  X  to  XII  show  the  variation  of 
probabilities  of  types  of  intervals  with  a  change  in  p, 
keeping  a  =  0.02. 


In  Tables  XIII  to  XV,  the  variation  of  the  probabilities 
of  types  of  intervals  are  shown  at  a  standard  deviation  (o) 
of  3.0,  whereas  in  Tables  XVI  to  XVIII,  the  value  of  standard 
deviation  is  altered  to  5.0.  Tables  XIX  to  XXVII  depict  the 
variation  of  the  types  of  intervals  with  change  in  sample 
sizes  at  different  values  of  p  and  a,  keeping  a  at  0.05. 
Finally,  Tables  XXVIII  to  XXX  show  the  variation  of  proba¬ 
bility  of  types  of  intervals  with  change  in  sigma  at  different 
values  of  p  and  n. 

In  all  these  tables  the  first  column  indicates  the  value 
of  the  parameter.  Columns  2,  3  and  4  give  the  probabilities 
of  occurrence  of  types  of  intervals.  Columns  5  and  6  show  the 
probability  that  the  predicted  value  is  contained  in  the 
interval.  Since  a  type  3  interval  is  of  the  form  (-«,+») 
the  probability  of  *2k+l  to  conta:*-ned  in  the  interval  is 
one.  Finally  column  7  verifies  the  relation  (3.23).  The 
attached  graphs  on  pages  58  to  62  show  the  trend  in  the  pro¬ 
bability  of  occurence  of  the  three  types  of  intervals  with 
changes  in  parameters  like  p,  a,  n  and  a. 

The  outcome  of  the  simulation  confirmed  the  following 
theoretical  conclusions : 

i)  The  probability  of  type  1  interval  increases  with  an 
increase  in  correlation  coefficient  p,  all  other  parameters 
remaining  constant. 

ii)  The  probability  of  type  1  interval  increases  with  an 
increase  in  sample  size,  keeping  all  other  parameters  constant. 
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iii)  The  probabilities  of  types  2  and  3  decrease  with 
an  increase  in  correlation  coefficient  or  sample  size. 

iv)  The  probabilities  for  the  four  types  of  intervals 
do  not  change  significantly  with  change  in  sigma. 

v)  The  probability  of  a  type  1  interval  increases  with 
an  increase  in  level  of  significance  a,  all  other  parameters 
remaining  constant. 

Section  6  deals  with  the  application  to  linear  trend 
model . 
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VI.  APPLICATION  TO  LINEAR  TREND  MODEL 

In  this  chapter,  the  results  of  chapter  3  are  extended 
to  stochastic  processes  that  follow  a  linear  trend  model. 
Let  {Xn}  be  a  discrete  stochastic  process  satisfying 


Xn  =  yn  +  en  * 


where  u 


n 


a  +  0n  , 


(6-1) 


(6-2) 


and  en  is  a  stationary  process  of  normal  variates  that 
satisfy  the  following: 


E(cn>  -  0 

f 

Var{e  }  =  a4 
n 


2  k 


Cor{en,£n+k)  =  a  p  ,  p  <  1 


(6-3) 


Such  a  process  is  called  a  model  for  linear  trend.  Krishnaih 
and  Murthy  [7]  showed  that  for  the  above  model  the  conditional 
distributions  of  X2k  given  X2k_-^  anc*  X2k+1  are  i-ndePendent 
normal  with 


E  ^X2k  I  X2k-1 ’  X2k+1^  "  S1  +  B2K  +  S3Xi 


and  Var  [X2k  |  X2k+]_]  =  a2(l-p2)  /  (1+p2)  , 


(6-4) 


(6-5) 


where 


*k  '  <X2k-l  +  X2k+1>/2 


1  .  -on  i.i.U. 
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In  section  A,  a  prediction  interval  for  X2k+^  is 
obtained  knowing  the  values  up  to  X2k  for  the  above  model. 

A.  PREDICTION  INTERVAL  FOR  X2k+1 

Using  regression  analysis,  a  prediction  interval  for 
X2k+1  can  be  obtained  tor  the  linear  trend  model,  knowing 
the  values  up  to  X2k>  Equation  (6-4)  can  be  written  as 

X2k  =  01  +  e2k  +  B3Xk  +  £k’  k  =  1'2*  • • • >  n, 
or  in  matrix  notation. 


X  2k  =  X  B  +  e  ,  e  v  N(0  ,  o2I) 

nxl  nx3  3x1  nxl 


The  design  matrix  X  is  of  the  form: 

nx3 


(6-6) 


1  1  (X1+X3)/2 
1  2  (X3+X5)/2 
1  3  (X5+X7)/2 


1  n  <X2n-l+!W'2 


and  B  is  the  column  vector  (3,,  30,  s»)T 
3x1  1  z  J 


The  estimator  B  can  be  obtained  from  the  normal  equations 
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rp  /v 

X  X  B 
3xn  nx3  3x1 


xT  X 

3xn  nxl^ 


or 


B  = 


-1  T 

*  —  — 2k 


where  S  = 


n 

In 

IX' 

n 

=  XT 

x  - 

in 

In2 

InX 

3xn 

nx3 

IX' 

n 

InXA 

i 

ix’: 

n 

(6-7) 


Since  S  can  be  shown  to  be  nonsingular  and  of  rank  3,  S~^ 
exists  and 


B  ^  N  (B,  a2S-1) 


(6-8) 


For  a  given  value  of  X', 

1x3  k 


=  [Ik  x^  ]  , 


the  predicted  value  X  =  X',  B 

lxlZk  1x3  k  3x1 


and  X2k  v  N  {X^  B  ,  a2(X^  )S'1(X^  )T} 


It  can  be  shown  that  the  residual 


x2k  “  x2k  ^  (1+s  ) }  > 


(6-9) 


where 


s2  =  X'  S_1(  X’  )T 
lxl  1x3k  3x3  3x1k 


As  an  example,  let  X  =  (x^ ,  x2 ,  . . . ,xlQ) 


Here  K  =  5  and  n  =  4. 


The  vector  X  01  =  (x0,  x, ,  x, ,  xe) 

4xlZk  2  4  6  8 
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and  the  design  matrix  X  = 

4x3 


(x-l+x3)/2 
(x3+x5) /2 

(Xj^+Xy) / 2 
(Xy+Xg) / 2 


X£  =[15  (x9+x11)/2]  . 

The  value  of  x-^  is  unknown.  A  prediction  interval  for  x^ 
can  be  developed  as  shown  below. 

From  equation  (6-9) 


t2k~  „  N  (o,  i) 


o/  1  +  s' 


(6-10) 


X2k  "  X2k 


1  +  s' 


=  t 


(6-11) 


The  quantity  t  has  a  student's  t  distribution  with  (n-3) 

2 

degrees  of  freedom,  where  R  is  the  usual  unbiased  estimator 
of  a2. 

At  a  given  level  of  significance  a  , 


(x2k  -  x2k)  -CT 

r_  n~+? 


=  1  -  a 


(6-12) 


The  term  in  the  brackets  of  the  left  hand  side  of  equa¬ 
tion  (6-12)  can  be  written  as, 
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(6-13) 


X2k  ‘  X2k 


4  Va/2 


and 


1  +  s  = 


i  +  [i  k  [s-1] [i  k  x-]T. 


The  expression  on  the  right  hand  side  can  be  simplified  to 
the  form 

y o  +  Ylxk  +  Y2xk2  • 

The  coefficients  yq>  Y-j_  and  y2  can  be  computed  since  k  and 
S  ^  are  known. 


Squaring  both  sides  of  equation  (6-13)  and  substituting  for 
2 

x?,  and  s  ,  it  can  be  simplified  to, 


2  2 

/N  /S  /N  t  ^  R 

('x2k~SI82k'03xk)  4  n-3  (yo+y lxk+Y2xk2) 


or 


t2  R2  t2  R2y 

(3  ^  2  °)  x,'2  +  (2(x  -8  -3  k)8  -  ct/ 2  o  1 »  i 

^3  y2  n-3  ;  *k  +  UCx2k  *1  B2  ;B3  n-3  ’  *k 

2  2 

A  /n  ry  t  /  r\  R  Y 

+  ((x2k-81-62k)2-  ^23°--  )  <  0  . 


above  inequality  reduces  to 


A  x’  +  Bx,'  +  C  <  0 


(6-14) 


where 


63  '  Y2Ca/2Ro/(n'3)  ' 

.  .  *  t2,7R2y, 

2(x2k'Bl'S2k)S3  '  a/n-3 — 
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(6-15) 


C  ^2k~  ^1~ 


2  2 
t  ,0r  r 

a/2  o  o 

~rn 


The  quadratic  equation  (6-14)  can  be  solved  for  x^Q  using 
the  formula. 


-B±/b2-4AC 
- 2A - 


if  A  ^  0  . 


Since  x^  =  %  (x2ic-l+x2k+l^  ’  tbe  Pre^iction  interval  for 
x2k+l  can  be  as 

D  -  E  <  x2k+1  <  D  +  E  (6-16) 

where 

D  =  ~x2k-l  '  B/A 

E  =  /(B2-4AC)/A  (6-17) 

2 

Note:  If  A=0,  then  no  interval  exists;  also  when  B  <  4AC 

E  becomes  imaginary  and  the  interval  is  taken  to  be  (-»,+«>). 

In  section  A,  a  numerical  example  is  given  to  verify  the 
model  developed  here,  with  the  help  of  simulated  samples. 


B.  NUMERICAL  EXAMPLE 

In  order  to  compute  the  prediction  intervals  for  a  future 
value  for  samples  satisfying  a  linear  trend  model  an  APL 
program  LINTREND  has  been  developed.  The  listing  of  this 
program  is  given  on  page  67.  Error  terms  having  the  correla¬ 
tion  structure  (6-3)  are  given  as  input  and  for  given  values 
of  a  and  8,  the  program  develops  the  vectors  X9k  and  X/. 
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For  a  sample  of  size  18,  the  output  of  the  above  program  is 
given  below. 

a  =  1.0  ,  8  =  0.5  ,  p  =  0.2  ,  N  =  18  ,  n  =  8. 

X  '  =  (2.094,  3.37,  4.056,  3.7595,  5.984,  7.848, 

8xlK  t 

7.821,  8,38s)1 


X  9,  =  (2.94,  3.983,  4.913,  2.685,  5.376,  7.849,  9.017, 
8xlzlc  T 

8 . 416) 1 


1 

1 


The  design  matrix  X 

8x3 


1 

1 

1 

1 

1 

1 


1 

2 

3 

4 

5 

6 

7 

8 


2.094 

3.37 

4.056 

3.759 

5.984 

7.848 

7.821 

8.388 


For  the  model  (6-4), 


B=(-0 . 147 ,  -0.597,  1.566) 
3x1 


The  coefficients  of  the  quadratic  equation  (6-14) 


A(X^)2  +  B (X£)  +  C  =  0  , 


are  given  by 


i-3 


A  =  1.44 
B  =  35.801 
C  =  -72.382 

and  the  values  of  D  and  E  are 
D  =  15.955 
E  =  28.561 

from  the  equation  (6-17)  , 
the  lower  limit  =  -  12.606 

the  upper  limit  =  44.516 

The  true  value  of  X-^g  (10.515)  is  contained  in  the 
interval . 
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I.  VARIATION  OF  TYPES  CF  INTERVALS  WITH  ROW 

SAMPLE  SIZE  =14 

SIGMA  =  1.3 

LEVEL  OF  SIGNIFICANCE  =  0.05 


ROW 

TYPE  1 

TYPE  2 

TYPE  3 

PRCB  1 

— - 

PROS  2 

CCN  REG 

0.100 

0.193 

0  .066 

0.741 

0.917 

C.  576 

0.  556 

0.20C 

C.342 

0.080 

0.  578 

0.  515 

C.  775 

0.  553 

0.  50 C 

C.  269 

0.  C98 

C.633 

0.  883 

0.736 

0  .949 

0.  70  0 

0.  369 

0.119 

0.512 

0.902 

0  .374 

0.549 

0.900 

0.539 

0 .128 

0.333 

0.939 

0.930 

0.  556 

11.  VARIATION  OF  T Y D ES 

OF  INTERVALS  WITH  ROW 

SAMPLE  SIZE 

=  22 

SIGMA 

=  1.0 

LEVEl  CF  SIGNIFICANCE  =  0. 

05 

ROW 

*YPE  1 

TYPE  2 

TY °E  3 

RR08  1 

RR9B  2 

CON  REG 

0  .100 

C.492 

0.  048 

C.  46  C 

0.  935 

C.438 

0.541 

C.  30  C 

C.  516 

0.  C77 

0.407 

0.  924 

0.714 

0.939 

0.  500 

0.483 

0.  110 

0.407 

0  .909 

0.818 

0.936 

0.700 

0.791 

0  .071 

0.138 

0.  934 

0.873 

0.  535 

0.900 

C.  902 

0.C50 

0.048 

0.  9<t7 

C.  920 

0. 54  8 

CON  REG  «  CONFIDENCE  REGION 
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III.  VARIATION  OF  TY  P  ES  CF  INTERVALS  W  IT  H  ROW 

SAMPLE  SIZE  =  30 

SIGMA  =  I. 3 

LEVEL  OF  SIGNIFICANCE  *  0.05 


ROW 

TYPE  1 

TYPE  2 

TYPE  3 

PROB  1 

FRQB  2 

CCN  REG 

0.100 

0.554 

0.036 

0.410 

0.960 

0.  722 

C.  56  8 

0.30C 

C.  861 

0.  033 

C.  106 

0.  940 

C.  758 

C.  54  0 

0.500 

C.  874 

0.  C44 

0.082 

0.  952 

0.341 

0.951 

0.700 

0.914 

0.  028 

0.058 

0.951 

0.9  29 

0.553 

0.900 

0.979 

0  .003 

0.013 

0.  952 

1.000 

0.  553 

IV. 

SAMPLE 

SIGMA 

VARIATION  OF  T 

SIZE  =  54 

=  1.  0 

YPES  OF  INTERVALS 

WITH 

ROW 

LEVEL  OF  SIGNIFICANCE  =  0.05 


ROW 

TYPE  1 

r  Y  PE  2 

mm 

PROS  1 

PROB  2 

CON  =EG 

0.100 

C.5  85 

0.  041 

0.  374 

0.  955 

C.  610 

C.  56  0 

0.30C 

C.  68  3 

0.  C66 

0.251 

0.  943 

0.773 

0 .946 

0.500 

0.718 

0.077 

0.205 

0.937 

0.766 

0.537 

0.700 

0.950 

0  .013 

0.037 

0.942 

0.  346 

0.  543 

0.900 

C.  992 

0.004 

0.  004 

0.  948 

1.  000 

C.  548 
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V.  VARIATION  0=  rypES  OF  INTERVALS  WITH  ROW 


SAMPLE  SIZE 

=  38 

S  I GMA 

=  1.0 

LEVEL  OF  SIGNIFICANCE  =  0.05 


ROW 

TY  PE  l 

T  Y  FE  2 

TY  P  E  3 

PROS  1 

PROB  2 

COM  REG 

0.100 

C.703 

0.  030 

0.  267 

0.  945 

C.  467 

0.  54  8 

0.30C 

C.  809 

0.  C36 

0.155 

0.  954 

0.806 

0  .956 

0.500 

0.847 

0.041 

0.112 

0.956 

0  .854 

0.957 

0.700 

0.966 

0  .015 

0.019 

0.959 

0.933 

C.  555 

0.900 

C.994 

0.C04 

0.002 

0.  95  7 

1.000 

C.  557 

VI 

.  VA3  IATION  OF  TYPES  OF 

INTERVALS  WITH 

ROW 

SAMPLE  SIZE 

=  1 2 

SI  GM  A 

=  1.0 

LEVEL 

CF  SIGNIFICANCE 

=  0.05 

ROW 

type  1 

TYPE  2 

PROS  1 

PROB  2 

CON  REG 

G.  IOC 

C.  573 

0.  C47 

C.330 

0.932 

0.723 

0.948 

0.300 

0.719 

0.  044 

0.23  7 

0.947 

0.773 

0.552 

0.500 

0.853 

0  .047 

0.10  0 

0.943 

0.  894 

0.  546 

0.70C 

C.980 

0.007 

0.0  13 

0.  543 

C.  857 

C.  54  8 

C.90C 

C.  998 

0.  COi 

0.  001 

0.  94  7 

1.000 

0  .547 

. - .  ■  „  . 

.  _ 
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VII.  V4RI4TI0N  OF  TYPES  OF  INTERVALS  WITH  ROW 


SAMPLE  SIZE  =  1  4 

SIGMA  =  1.0 

LEVEL  CF  SIGNIFICANCE  =  O.IO 


ROW 

US 

TYPE  2 

TY°  E  3 

PROB  1 

PROB  2 

CON  REG 

C.10C 

C.  269 

0.  Ill 

C.  62  C 

0.  87  0 

0.468 

0.906 

0.300 

0.441 

0.  108 

0  .451 

0.346 

0  .667 

0 .89  6 

0.500 

C.380 

0 .138 

0.482 

0.858 

C.  732 

0.  SC 9 

0.700 

0.501 

0. 138 

0.361 

0.  85  8 

C.  790 

c.  soo 

C.  90C 

C.  666 

0.  114 

C.22  C 

0.  883 

0.868 

0.907 

VII  I  .  VARIATI  ON  0=  TYPES  CF  INTERVALS  WITH  RCW 

SAMPLE  SIZE  =22 

SI  GVA  =  1  .0 

LEVEL  OF  SIGNIFICANCE  =  0.10 


ROW 

TYPE  1 

TYPE  2 

TY°E  3 

PR  03  1 

FRC8  2 

CON  REG 

0.  100 

0.  524 

0.077 

0.399 

0.876 

0.468 

0.894 

0.300 

0.599 

0 .101 

0.300 

0.873 

0.  564 

0.  €8C 

0.50  0 

0.617 

0. 122 

0.261 

0.  861 

C.754 

0.884 

0.  70  C 

C.853 

0.  C66 

0.081 

0.  885 

0.773 

0.887 

0.900 

0.941 

0.  023 

0.031 

0.391 

0  .786 

0  .89  1 

IX.  VARIATION  0=  TYPES  CF  I  NTfcRV  AlS  WITH  ROW 


SAMPLE  SIZE  =  30 

SIGMA  =  I  .0 

LEVEL  OF  SIGNIFICANCE  =  0.10 


ROW 

TYPE  1 

TYPE  2 

TYPE  3 

PROB  1 

FRQB  2 

CON  REG 

0.  100 

0.595 

C.070 

0.335 

0.899 

0.571 

0.910 

0.300 

0.906 

0  .025 

0.069 

0.891 

0.  560 

0.  ESC 

0.500 

C.933 

0.C24 

0.043 

0.  902 

C.  875 

C.906 

C.70C 

C.  951 

0.  C22 

0.027 

0.902 

0.909 

0  .905 

0.900 

0.989 

0.006 

0.005 

0.897 

1  .000 

0.898 

X. 

VARI ATI CN  OF 

TY  P ES  OF 

INT  ERV  ALS 

W  ITH 

ROW 

SAMPLE 

SIZE 

=  L  4 

SIGMA  =  I  .0 

LEVEL  OF  SIGNIFICANCE  =  0.02 


ROW 

TYPE  1 

TYPE  2 

TYPE  3 

PR  CB  1 

FROB  2 

CON  REG 

0.100 

0.160 

0  .028 

0.812 

0.  969 

0.  571 

0.  983 

0.300 

0.254 

0.C42 

0.704 

0.  949 

C.738 

C.  976 

C.  50  C 

C.  169 

0.  C63 

C.  76  8 

0.929 

0.825 

0 .977 

0.700 

0.209 

0.088 

0.703 

0.933 

0  .932 

0.980 

0.900 

0.3  81 

0.121 

0.498 

0.961 

C.  975 

C.  982 

XI.  VARIATION  CF  TYPES  QF  INTERVALS  WITH  ROW 

SAMPLE  SIZE  =  22 

SIGMA  =  1  .a 

LEVEL  OF  SIGNIFICANCE  =  0.02 


ROW 

TYPE  1 

TYPE  2 

TYPE  3 

PR  CB  1 

FROB  2 

CCN  REG 

0.100 

0.464 

0  .026 

0.510 

0.  970 

0.  731 

0.  97  9 

0.300 

C.427 

0.054 

0.  519 

0.  96  7 

C.  870 

0. 979 

0.  50 C 

C.  338 

0.  C76 

0.586 

0.  938 

0.987 

0  .978 

0.700 

0.691 

0.  073 

0.2  36 

0.965 

0  .973 

0.974 

0.900 

0.82  5 

0.073 

0.10  2 

0.979 

0.986 

0.  982 

XII. 

VARIATION  HF  TYRES  OF 

INTERVALS  WITH 

ROW 

SAMPLE 

SIZE 

o 

rA 

II 

S  IG«* 

=  1.0 

LEVEL 

OF  SIGNIFICANCE 

=  0.02 

ROW 

"rY  PE  1 

mm 

TYPE  3 

PROS  1 

PRO B  2 

CON  REG 

0.100 

0.525 

0.  C19 

0.  456 

0.  983 

C.  737 

C.  986 

C.30C 

C.  830 

0.  C22 

0.148 

0.  972 

0.773 

0  .972 

0.500 

0.810 

0.  032 

0.158 

0.978 

0.938 

0.98  0 

0.700 

0.846 

0  .051 

0.10  3 

0.979 

C.961 

C.  98C 

0.90C 

C.961 

0.015 

0.024 

0.  977 

1.000 

0.  97  8 
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X  1 1 1  .V  AR  I  AT  ION  n-  ^YPES  OF  INTERVALS  WITH  ROW 

SAMPLE  SIZE  =  12 

S  IG*A  =3.0 

LEVEL  OF  SIGNIFICANCE  =  0.05 


ROW 

TYPE  l 

TYPE  2 

TY°E  3 

PROS  1 

PROS  2 

CON  REG 

0.100 

0.06  7 

0.  C78 

C.855 

0.  851 

C.  500 

0.951 

C.  30  C 

C.  179 

0.  125 

C.  696 

0.  855 

0.704 

0.937 

0.500 

0.445 

C.  117 

0.438 

0.899 

0  .829 

0.935 

0.700 

0.736 

0.083 

0.176 

0.938 

C.  898 

C.  545 

o.yoc 

C.901 

0.  051 

0.048 

0.  948 

C.  941 

C.550 

XIV.  VARIATION  OF  T 

YPES  OF 

I NTERVALS 

WITH 

ROW 

SAMPLE  S  ize 

=  26 

sigma 

=  3.0 

LEVEL  CF  SIGNIFICANCE 

=  0.05 

ROW  TY°E  1 

TYPE  2 

type  3 

PROS  1 

PROB  2 

CON  REG 

C.  IOC  C.  652 

0.  C43 

0.305 

0.  936 

0.651 

0.943 

0.300  0.255 

0.106 

0.639 

0.929 

0.699 

C.550 

0.500  0.561 

0 .110 

0.329 

0.933 

0.855 

C.  545 

0.70C  C.841 

0.  062 

0.09  7 

0.  941 

C.  903 

0.  544 

0.90C  C.  962 

0.  C16 

0.022 

0.  946 

0.938 

0.947 

— - —  —  —  _  _ 

—  _  — —  _ 

....  ————  —  — 

XV. VAR  IATION  0:‘  TYPES  OF  INTERVALS  WI^H  ROW 


SAMPLE  S  IZE  =  33 

SIGMA  =3.0 

LEVEL  CF  SIGNIFICANCE  =  0.05 


ROW 

TYPE  1 

TYPE  2 

T  Y°  E  3 

PROS  1 

PROS  2 

COM  REG 

C.  IOC 

C.  075 

C.  C71 

C.854 

0.  88  0 

0.620 

0  .964 

0.300 

0.374 

0.  Ill 

0.515 

0.922 

0  .721 

0.940 

0.500 

0.674 

0  .091 

0.235 

0.942 

C.  857 

C.  946 

0.70C 

C.  8S7 

0.034 

C.  069 

0.  95  C 

C.  912 

0.  952 

0.90C 

C.  979 

0.  C07 

C.  014 

0.  952 

1.000 

0  .953 

XVI  .VARI  ATI  ON  OF  TYPES  CF  INTERVALS  WI^H  OQW 

SAMPLE  SIZE  =  22 

SIGI*A  =5.0 

LEVEL  CF  SIGNIFICANCE  =  0.05 


ROW 

TY°  E  1 

TY»E  2 

TY°E  3 

PROS  1 

PR3B  2 

CCN  REG 

0.  100 

0.067 

0.078 

0.855 

0.851 

0  .500 

0.951 

0.300 

0.182 

0  .127 

0.69  1 

0.863 

C.717 

C.  939 

0.50  0 

C.  444 

0.  115 

0.  441 

0.  901 

C.  317 

0.  93  5 

0.  70 C 

C.  716 

0.  102 

0.182 

0.  936 

0.892 

0  .943 

0.  900 

0.901 

0.052 

0.047 

0.948 

0.942 

0.950 
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XVII  .VARI  ATI  ON  0=  TYPES  CF  INTERVALS  WITH  ROW 


SAMPLE  SIZE  =26 

SI  GMA  =5.0 

LEVEL  OF  SIGNIFICANCE  =  0.05 


ROW 

TYPE  i 

TYPE  2 

TYPE  3 

PROS  1 

PROS  2 

CON  REG 

0.  100 

0.193 

0.097 

0.710 

0.896 

0  .619 

0.943 

0  .30  0 

0.255 

0.107 

0.638 

0.929 

C.701 

C.  95C 

0.5  OC 

C.56  0 

0.  Ill 

0.  329 

0.  938 

C.  856 

0. 949 

0.  70 C 

C.  841 

0.  C62 

C.097 

0.941 

0.903 

0  .944 

0.  900 

0.962 

0.016 

0.022 

0.9*7 

0.938 

0.948 

XVIII 

. VARI A 

T I  ON  OF  type-5  if 

IMT  ERVALS 

WITH 

ROW 

SAMPLE 

SIZE 

=  30 

SIGMA 

=  5  .0 

LEVEL 

OF  SI3 

NIPICANCE 

=  0.  05 

ROW 

TYPE  1 

TYPE  2 

TYPE  3 

PRCB  1 

c  PD  3  2 

CCN  REG 

0.100 

0.075 

0  .075 

0.850 

0.89  3 

C.  640 

C.  965 

0.300 

C.350 

0.  110 

C.  540 

0.  92  0 

C.  709 

C.  940 

0.  50 C 

C.  664 

0.  C97 

C.239 

0.  950 

0.845 

0.952 

0.700 

0.897 

0.034 

0.069 

0.950 

0.912 

0.952 

0.900 

0.979 

0  .007 

0.014 

0.  952 

1.000 

0.  953 

— .  — —  — - 

.  —  _ _ _  _  —  —  _ 
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XIX.  VARIATION  OF  ”YPES  CF  INTERVALS  WITH  N 


CORR.COEFF.  =  0.300 

SIGMA  =  1.0 

LEVEL  OF  SIGNIFICANCE  =  0.05 


N 

TYPE  1 

TYPE  2 

TYPE  3 

PRCB  1 

FROB  2 

CCN  REG 

4 

G  .3  00 

0  ^040 

0.660 

0.  867 

C.600 

0.  544 

6 

C  .342 

0.  C80 

0.  578 

0.  515 

C.775 

0.  553 

8 

C.  587 

0.  C57 

C.356 

0.  942 

C  .73 7 

0.951 

10 

0.516 

0.077 

0.407 

0.924 

0  .714 

0.539 

12 

0.259 

0 .106 

0  .635 

0.  946 

C.  698 

0.  554 

14 

C.861 

0.C33 

0.  106 

0.  940 

C.  758 

C.  54  C 

XX 

.  VAR  IATION  OF  TYPES  OF 

I NTERVAL 

S  WITH 

N 

C3RR. 

COEFF. 

=  0  .500 

SIGMA 

=  1.0 

LEVEL 

CF  SIGNIFICANCE 

=  0.05 

N 

ryoE  i 

TYPE  2 

TYPE  3 

PR  08  1 

PROS  2 

CON  REG 

4 

C.232 

0.  C49 

0.719 

0.  353 

0.694 

0.95  1 

6 

0.269 

0.098 

0.633 

0.883 

0 .786 

0.549 

8 

0.540 

0  .083 

0.377 

0.  935 

C.319 

C.  55C 

10 

C.*83 

0.  110 

C.407 

0.  909 

C.  818 

C.  536 

12 

C.  561 

0.  109 

C.33C 

0.  941 

0.353 

0.55  1 

14 

0.874 

0.  04^ 

0.082 

0.952 

0  .841 

0.551 
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XXI  .  VAP.I  ATI  ON  0=  TYPES  CF  INTERVALS  MTn  M 


CQRR  .COEFF .  =  0. 70C 

SIGMA  =  1  .0 

LEVEL  CF  SIGNIFICANCE  =  0.05 


N 

TYPE  1 

TYPE  2 

TYPE  3 

PRQB  1 

PR3B  2 

CON  REG 

4 

0.308 

0.07? 

0.620 

0 .39o 

0  .861 

0.958 

6 

0.36  9 

0 .119 

0.512 

0.  9 02 

C.  874 

C.  948 

6 

C.586 

0.  095 

C.  319 

C.  920 

C.  884 

C.  842 

10 

C.  791 

0.  C71 

C.  138 

0.  934 

C  .873 

0  .939 

12 

0.840 

0.065 

0.095 

0.940 

0  .908 

0.844 

14 

C.914 

0  .028 

0.058 

0.951 

C.9  28 

C.  953 

XXII  .VAR  I  AT  ION  OF  r 

Y  °  ES  CF 

INTERVALS 

WITH 

N 

CORR. 

CCEFF. 

=  0.300 

S  IGMA 

=  3.0 

LEVEL 

CF  SI  G 

NI  FI  CANCE 

=  0.05 

N 

"YPE  L 

TYPE  2 

TY°E  3 

=>ROB  1 

PRQB  2 

CCN  REG 

4 

C  •  Go  7 

0.  050 

0.  383 

0.  642 

C.  580 

C.  855 

6 

C.  090 

0.  C93 

C.  317 

0.689 

C.710 

0.845 

8 

0.262 

0.  085 

0.653 

0.924 

0  .718 

0.856 

10 

0.179 

0 .125 

0.696 

0.855 

C.  704 

C.  937 

12 

C.255 

0.  106 

C.  638 

0.  929 

C.  698 

C.  95  C 

14 

C.  374 

0.  Ill 

C.515 

0.922 

0.721 

0.84  0 

rr-M  —  -  - 

— _ ■  ->  _  _  . 

53 


XXIII  .VARIA  TION  OF  TYPES  OF  INTERVALS  WITH  N 

CORR.CCEFF.  =  0.500 

SIGMA  =3.0 

LEVEL  CF  SIGNIFICANCE  =  0.05 


N 

TYPE  1 

TYPE  2 

TYPE  3 

PROS  1 

PROS  2 

CON  REG 

4 

C.  108 

0.  C68 

0.824 

0.  713 

0.750 

0.95  2 

6 

0.226 

0.105 

0.669 

0.863 

0.781 

0.946 

3 

0.376 

0.113 

0.506 

0.923 

C.  822 

C.  95 C 

10 

C.445 

0.  117 

0.  438 

0.  699 

C.  829 

C.935 

12 

C.  561 

o.  no 

C.329 

0.  938 

0.355 

0.949 

14 

0.674 

0.091 

0.2  3  5 

0.942 

0.857 

0.948 

XXIV. 

.  VARIATION  OF  TYFES  OF 

INTERVALS  WITH 

M 

CORR.CQEFF. 

=  0.  70 C 

SIGNA 

=  3  .0 

LEVEL 

OF  SIGNIFICANCE 

=  0.  05 

N 

TYPE  1 

TYPE  2 

TYPE  3 

PROS  1 

PROB  2 

CCN  REG 

4 

0.216 

0  .077 

0.707 

0.829 

C.  870 

C.  952 

6 

C  .354 

0.  122 

0.  524 

0.  898 

C.  885 

C.  950 

8 

C.  557 

0.  118 

C.325 

0.  92  6 

0.864 

0.94  3 

10 

0.736 

0.088 

0.176 

0.938 

0  .898 

0.945 

12 

0.841 

0  .062 

0.097 

0.  941 

C.903 

0.  944 

14 

0.897 

0.  C34 

C.  069 

0.  950 

C.  912 

C.  952 

-r— r-_^l 

_  — n 
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XXV.  VARIATION  OF  TYPES  OF  INTERVALS  WITH  N 


CORR.CCEFF.  =  0.300 

S  I  GMA  =5.0 

LEVEL  OF  SIGNIFICANCE  =  0.05 


N 

TYPE  1 

TY  FE  2 

TY° E  3 

PRO0  1 

PROB  2 

CON  REG 

4 

C.  06  8 

0.C51 

0.  88  1 

0.  632 

C.  588 

0.  954 

6 

C.  088 

0.  C94 

C.  818 

0.  682 

0.713 

0.945 

8 

0.138 

0.  093 

0.769 

0.884 

0.710 

0.957 

10 

0  .182 

0 .127 

0.69  1 

0.86  3 

C.  717 

C.  939 

12 

C  .25  5 

0.  107 

C.  628 

0.  925 

C.  701 

C.  950 

14 

C.  350 

0.  110 

C.  540 

0.  92  0 

0.709 

0.940 

XXVI 

.  VARIATION  OF  TYPES  CF 

I NTERVAL  S 

WTh 

N 

CO  PR  . 

COEFF. 

=  0. 50  C 

SIGMA 

=  3.0 

LEV  EL 

CF  S 1GN1F  ICAVCE 

=  0.05 

N 

TYD  E  1 

TYPE  2 

TYPE  3 

PROB  1 

PROB  2 

CON  REG 

4 

0.107 

0.  C69 

0.824 

0.701 

0  .739 

0.950 

6 

0.215 

0 .109 

0.676 

0. 860 

C.  761 

C.  944 

8 

C.363 

0.  126 

0.  511 

0.  928 

C.  81  0 

C.  95  C 

10 

C.  444 

0.  115 

C.441 

0.  901 

0.817 

0.935 

12 

0.560 

0.111 

0.329 

0.938 

0  .856 

0.949 

14 

0.664 

0  .097 

0.239 

0.  950 

C.845 

C.  952 

.  _  __  _  _  _ 

..  ___ _ ___ _ 

_ r 

. - r-.  .  _  . 
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XX  VII  .VARI  ATI  GN  CF  TY  FES  OF  INTERVALS  WITH  N 


CORR.COEFF.  =  0.900 

SIGMA  =  5  .0 

LEVEL  OF  SIGNIFICANCE  =  0.05 


N 

TYPE  1 

TYPE  2 

TYPE  3 

PROB  1 

F ROB  2 

CON  REG 

4 

0.237 

0  .103 

0.66  C 

0.873 

C.922 

0.  562 

6 

0.53  7 

0.  129 

C.  334 

0.  93  7 

C.  930 

C.  557 

3 

C.  777 

C.  C80 

0.  143 

0.943 

0.888 

0.947 

10 

0.901 

0.  052 

0.047 

0.948 

C  .942 

0.550 

12 

C  .96  2 

0  .016 

0.022 

0.947 

C.  538 

C.  546 

14 

C.S79 

0.  007 

0.  014 

0.  552 

1.000 

0.  553 

XXVIII  .VARIA  TION  OF  INTERVALS  WITH  SIG^A 

SAVDLE  SIZE  = 

CORR.COEFF.  =  0.700 

LEVEL  CF  SIGNIFICANCE  =  0.05 


SIGMA 

ty°e  1 

TY  D  E  2 

T  Y0  E  3 

PROS  1 

PROB  2 

CON  o  EG 

1 

C.  369 

0.  119 

C.  512 

0.  902 

0.374 

0.54  9 

2 

0.352 

0.  123 

0.525 

0.901 

0  .886 

0.95  1 

3 

0.354 

0.122 

0.524 

0.898 

C.  885 

C.  550 

4 

C  .354 

0.  121 

0.  525 

0.  658 

C.  884 

C.  550 

5 

C.  353 

0.  122 

C.  52  5 

0.  398 

0.885 

0.950 

XXIX.  VARIATION  of  INTERVALS  with  sigma 


SAMPLE  SIZE  =  18 

CORR.CCEFF.  =  0.500 

LEVEL  QF  SIGNIFICANICE  =  0.05 


SIGMA 

TYPE  1 

TYPE  2 

TYPE  3 

PROS  1 

FROB  2 

CCN  REG 

1 

C.540 

0.083 

0.377 

0.935 

C.819 

C.  95  0 

2 

C.3  96 

0.  116 

0.  488 

0.  919 

C.  819 

0.  94  7 

3 

C.  376 

0.  118 

C.  5  06 

0.923 

0.322 

0.950 

4 

0.372 

0.118 

0.51  0 

0.922 

0.822 

0.950 

5 

0.363 

0 .126 

0  .511 

0.  928 

C.810 

C.  950 

XXX 

.VAR  I  A 

TION  3F 

INTERVALS 

WITH  SIGMA 

SAMPLE 

SIZE 

=22 

COR°  .me FF. 

=  0. 

500 

LEVEL 

QF  SIS 

NI  F  I  CANC  E  =  0.05 

SIGMA 

’’Y  PE  1 

mu 

TY  DE  3 

DR0  8  1 

RROB  2 

CON  REG 

1 

C.483 

a 

H 

H 

• 

o 

0.  40  7 

0.  909 

C.  818 

C.  936 

2 

C.  446 

0.  115 

C.439 

0.  901 

0.826 

0.936 

3 

0.445 

0.  117 

0.43  8 

0.899 

0.828 

0.935 

4 

0.446 

0.115 

0.439 

0.  899 

0.826 

0.  935 

5 

0.444 

0.  115 

0.441 

0.  901 

C.  817 

C.  935 

- r - r  — 

_  — —  — 

_  —  —  _ 
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SAMPLE  SIZE  =  30 

SIGMA  *  1.0 

SIGNIFICANCE  LEVEL  =  0.05 

1,  2  and  3  on  the  curves  indicate  the  types  of  Intervals 
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PROB.  OF  TYPES  OF  I 


TYPES  OF  INTERVALS  VS  N 


SIGNIFICANCE  LEVEL 


TYPES  OF  INTERVALS  VS 


NT. 


5.  PROB.  OF  TYPE  1  INT.  VS  N 
(VARIATION  WITH  a) 


SIGMA 


I 


C  PROGRAM  TO  CALCULATE  PROBABILITY  CF  TYPE  1,TYPE  2 ,TY  d  E 

C  3t T  YP  E  4  INTERVALS  FpR  SIMULATED  SAMPLES. 

C  PROGRAMMER  T  S  MURTHY  SEP  1979. 

C  *#*#  ###*!!:**#*#  **********  **£**:)!**:*#**$:*  *$**#$>■.$*:$>>:> 

Cl  MENS ION  Z(55),X(55  ),V(100,55),XI150),YI(  50) ,S<  55) , 
1IC  (5),  IP.  (  5),  IPC(  5)  ,  STA  T(  10,  7)  ,  VS  (  10 , 10 ,  7)  ,  I  VS  ( 1  C) 

CALL  OVFLOW 
INDEX  *1 
SIGMA=1.0 

1  READ! 5,2)  K,T 
WR ITE  ( 6 ,2 ) K, T 

2  FOR  *AT(  IX,  12,  2X  =F  5.  3) 

IF ( K  .EQ.  0  )  GO  TO  460 
ROW =0.10 
DO  300  13=1,  5 
DO  10  J=1  ,5 
8  =  ( 1-ROW*  *2  )**.5 

C  SIMULATION  OF  SAMPLES 

I SEED= 12345 

10  IPCIJ)=0 

CO  250  I  A  =  1 , 10 
DO  5 C  M=l ,  100 
CALL  SNORM(ISEED,Z,K) 

X  ( 1 )  =S  IGMA*  Z  ( 1 ) 

CO  30  J=2 , K 

30  X(J  )=ROW*X(  J- 1)*3  *S  IGMA*Z(  J) 

DO  40  L  =1  , K 

40  V(M,L)=X(L) 

50  CONTINUE 

DO  60  II =1,5 

60  I R  ( 1 1 )  =0 

CO  200  1=1,  IOC 
DO  70  J=1,K 
S(  J)  =V(  I  ,  J) 

K=  K  -5 
KX=K/2 
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70 


80 


00  EC  L-l.KK 
YI(U«S(2*U 
N=KK-1 

DO  90  LL=1,N 
90  XI  (  LL)=(  S(2*LL-1)  +  S(2*LL+1J  1/2.0 
XSUMO.O 
YSUM=0.0 

sxx=c.  o 

SXY=0.0 
SYY=C.0 
DO  100  KL  =  1»  5 
100  IC ( KL  J  =0 

DO  110  M=1,N 
YSUM=YSUM+YI ( Ml 
XSU  M=XSUM  *X  I  (  M 1 
110  CONTINUE 

X8  =  XSUM/N 
Y9=YSUM/N 
DO  120  M=1,N 
$XX=SXX+<  XI  (  M )  — X  8 )  **2 
SYY  =  SYY + ( Y I l M  1  -Y 3) **2 
SXY=SXY+(  XI ( V )-X8) *{ YI (  MJ-YB) 

120  CONTINUE 

VRE$=(SYY-{  (S  XY**2  l/SXX))/(N-2) 

EH= SXY/SXX 

AH  =  YB-BH* XB 

SS=(T**2  ) * V  R  ES 

A= ( EH**2  l-SS/SXX 

°=S<2*KK)-Y9 

0*-  (  2*9  H*  P ) 

C3 ( °**2 1- (S  S  * ( N+l  )  }/N 
F= ( B  H**2  )*SXX/SS 
$SQ31.0-N*<  P**2)  /(  SS*(N+1J) 

IF  (F  .LT  .  SS  C)  GO  -n  15Q 
IF  (  F  •  EO  .1  «0  )  GO  to  50  0 
D*(  2.*XB  l-S(K-l)  +2.  *9H*D/A 
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E=( 2. /Al* (( ( SS*(  P**2) /SXX)+(N+1)*A*SS/N)**.5) 

PI L=D-E 
P IR=0+E 
P  VA  L=S(  K+  1) 

IF  C  F  .GT.  1.0)  i  0  tq  130 
IC  (2  1=  IC  (  2  1  +1 

IF( PVAL.LE.PIR  .OR.  PVAL.GE.PIL)  IC( 5) =  IC (5 )+l 
GO  T0  160 

j  130  1C  (1  )=  IC <  1  )+l 

IF(OIL.LE.PVAL.AMD.PIR.GE.PVAL)  IC(4)=IC(4)+1 
P.  GO  TO  160 

>  150  JC  <3  1  =  IC(3  1 +1 

160  01  170  J=  1,  5 

170  IR(  J) =1  R( J) +1 C( J) 

K=K*5 

|  200  CONTINUE 

j  CO  220  J  =  lf  5 

220  IPC( J)=IPC( JJ+IRI JJ 
250  CONTINUE 

C  PRINT  STATIST  ICS 

ST AT(  IR,1)  =  P0U 
STAT(  IB  ,21=1  PCIU  /I  CCO.  C 
ST  A” ( I B , 3 1  =  1  PC (2 ) /I 000 . 0 
ST  AT  (  I B,  4  )=  IP  C  ( 3  J/1000.0 
IF(  STAT(  16,2)  .EQ.O.  C)  GO  TO  275 
STATt IB ,5)=IPC(4)  /< ST AT ( I B»  2 ) *  1 00 0 .0 ) 

GO  TC  280 

275  ST*T(IB,5  J=I°C<4) 

280  IF< STAT( IB,3)  .EQ.  0.0)  GO  TO  285 

STAT(I3,6)=IPC(5)  /(  ST  AT  (18,31*1000.0 
GC  TC  290 

285  ST  AT  ( 18,6  1=  IpC(5  1 

290  SI G=  STA T(  16  ,2  )  *S  TA  T  (  I9,5)  +  STAT(IB,31*STAT{  IB, 6) 
1  +STAT( IB  ,4) 

ST  AT<  18,7  1=S  IG 
R0W=R0W+0  .2 
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300  CONTINUE 

00  305  I Q=1 ,5 
CO  305  JQ  =  1 »  7 

VS ( INDEX, IQ  ,JQ)  =  STAT(  IQ,JQ) 

305  C JNTI NUE 

IVS(  INDEX )=N 
INDE  X=INDcX+ 1 
I  5  Z  =  K-5 

WRITE  (6, 3  10)  ISZ,  S  IGM4,N 
WR  ITE  (  6,325) 

MI  TE(  6, 3  SOM  (STAT<  K,U  ,L=1  ,7)  ,K  =  1  ,5) 

310  FORMAT  (IX  »'  SAMPLE  SIZE  =  ',15,'  SIGMA  =  *,F5.0,  ' 

IN  =  '  ,  15,  /) 

325  FORMAT!  IX,'  ROW  TYPE  1  TYFE  2  7Y  ^E  3 

1FRGE.1  FRCB.2  CGN. 0  EG  *,/,70('-')) 

350  FORMAT  (  7(  F8  .3,  2X  ),/  ) 

WR I TE ( 6 , 485 ) 

GO  TO  1 
460  CCR  F=0 . 1 

INCEX=  INDEX-1 
CO  470  J  =  1 , 5 

WRITE  (  6,472)  CORR, SIGMA 
WR  Ir  E  (6  ,475  ) 

DC  471  1=1,  INDEX 

WRI7E(e>»480)(IVS(I)  y(VS(I,J,Ki  »K  =  2  ,7 ) ) 

471  CON'INUE 

COR  F=  CORR  +0  .2 
470  CONTINUE 

472  FORMAT!  '  COPR.COEFF.  =',F 5.3,'  SIGMA  *  ',F5.0,//) 
475  FORMAT ( '  SAMPLE  SIZE  “ Y  P  E  1  TY  3  E  2  ty3E  3  °R  0 

15.1  °R  DR.2  COM. DEG  ' ,  / ,  70  (  •-  ' )  ) 

480  FORMATd  5,5X,6(F8.  3,2X)  , /) 

500  STOP 
ENC 
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APL  PROGRAM  LINTREND 


V  LINTREND  X 

Cl]  5«-l+(N+0,5X  (  \NffX)  ) 

C2]  ;:f  <-::cn-i  j 

C3]  x<.(N-2)tx 

C4]  K«-<f*>x0.5 

C5]  ><i<-lO 

C6]  «<-l 

C7]  ->11X  \  K  <M«-M+1 

C8]  5<1<-<«C(  <  2xM )  -3 )  ]+xC  (  <2XM)-1)])X0.5 

C9] 

CIO]  -»7 

Cll]  fi<-xc2x\K] 

C 12]  N«-K-l 
C 13]  'rifWfn 

c 1 4 ]  1  THE  VECTOR  XI  IS  *JXI 
CIS]  '  THE  vector  TI  IS  ' }YI 
C16]  '  the  value  of  xf  is  • ;xp 

C 1 7]  cn-Nfi 
C 18]  C2<-IN 
C  1 9  ]  C3<-XI 

C20]  DMT<-(3f  N)f  Cl,  C2,C3 
C21]  DM«-$DMT 
C22]  SfDMT+.xDM 
C23]  sifgs 

C24]  8Hf(5I+,XDMT)+,X'n 

C25]  SHSf(+/( (TI ) *2) )-( (HTI+, XDM)+, XBH) 

C26]  A04-SIC1  >  i:  +  (N+i )  X  <  (2XSIH1  )2]  )  +  (><  + 1 )  XSIC2J2]) 
C27]  «0<-AO+l 

C28]  Ai<-<2xsicir33)+2xSi[;2»3]x<N+i ) 

C29]  «2*-siC3f3] 

C30]  CFf4XSHSr(N-3) 

C31  ]  A4-(bhc3]*2)-<02xcf) 

C32]  E««-(2xeHt;3]xwv(XC2x  i<]-bhci]-bh|;2]x  k))-aixcf 
C333  c*.<  w«2)-«0xcf 

C34]  4iMxi(8x2U4xflxc 
C35]  £♦-<  (d*2)-4xaxc)  *0.5 

C  36]  e<-eh.a 

C37]  D4--1XXC  <2XK  )-13-B-rA 
C38]  1  THE  LEFT  LIMIT  IS  ' fD-E 
C39]  '  THE  RIGHT  LIMIT  IS  ' }D+E 
C40]  -»0 

C413  IM* '  NO  interval  exists  ' 

V 
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APL  PROGRAM  PREDICT 


Cl]  n  y'  15  ASSUMED  TO  BE  A  VECTOR:  OF  EVEN  NO  OF  ENTRIES 

C23  K«-(f»)x0.5 

C33 

C4]  *<-l 

C53  ->9XlK<MfM+l 

C6 ]  *1<-<KC  <  (2xm>-3)  ]  +  ::[  <  < 2xM ) -1 ) 3 ) xO . 5 

C7]  xifxi,xi 
CS3  -»5 

C9]  Tl4-i<C2X\K] 

CIO]  h<-k-1 
C  1  1  3  TlfNfTI 
C 123  TB<.<+/TI)iM 
C  133  XB«-<+/xx  )tN 

C 1 43  SXff+/(  (XI-XB)x( YI-TB)  ) 

C153  SKXf+/( (XI-XB) *2) 

C  1 63  ST T<-+/(  (YI-YB)  *2) 

C173  BHfSXY+SXX 
C183  AHf.YB-BHxXB 

C  1 93  VRESf  (SYY+cix  (  (s:ct  *2)tS::::)  )  )+n~2 

C203  g<-  12.706  4.303  3.182  2.776  2.571  2.447  2.365  2.306 

C213  g<-g»  2.201  2.179  2.16  2.145  2.131  2.12  2.11  2.101  2.09 

C223  g<-g»  2.08  2.074  2.069  2.064  2.06  2.056  2.052  2.048  2 

C233  g«-g,  2.041  2.04  2.04  2.04  2.035  2.03  2.03  2.03 

C243  CfG,  2.025  2.025  2.02  2.02  2.02  2.02  2.015  2.01  2.01 

C253  t«-gcm-23 

C263  fl<-(BH*2)-(  (rSXX)  XSf.(  (T*2)  XVRES)  ) 

C273  B<-~lX(2xBHxF<-(  >SC2X  K  3-  TB  )  ) 

C283  ckp<2)-5x(('<+1)tN) 

C293  F«-<BH*2)  XSXXiS 

C303  1  THE  F  STATISTIC  IS  1  }F 

C31 3  ssq<-i-(Ht<n+i  )  )  x<p*2>+s 

C323  ->T1X\F>1 

C333  -»T2x  \ f=i 

C343  -»T3x  \ f < ssq 

C353  1  THE  INTERVAL.  IS  OF  TYPE  2  ’ 

C363  A0 
C373  TitA\0 

C383  »«-<2xn»>-xe  (2xK)-i]+2xbhxp-tA 

C393  Ef(rA)x2X(5x( ( ( F*2 ) tSXX )  +  < (M+l)rH)xA) )*0.5 

C403  1  the  left  limit  is  1  ;d-e 
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